Interaction of 3-level atom with radiation 



o 
o 

< 



(N 



L.Accardi, K.Imafuku, S.V.Kozyrev 
February 1, 2008 



Abstract 



The interaction of 3-level system with a quantum field in a non-equilibrium state is 
considered. We describe a class of states of the quantum field for wich a stationary state 
drives the system to inverse populated state. We find that the quotient of the population 
of the energy levels in the simplest case is described by the double Einstein formula which 
involves products of two Einstein emission/absorption relations. Emission and absorption of 
. radiation by 3-level atom in non-equlibrium stationary state is described. 

3t ! 1 Introduction 
o 

In the present paper we consider a 3-level system (for example an atom) interacting with radiation. 
In the work [JTJ] we apply the technique developed here to describe stimulated emission for 3-level 
O-i! atom interacting with radiation in a non-equilibrium state. In the work [[| a 2-level quantum 
system interacting with quantum field was considered. In a stationary nonequilibrium state 
for an n-level system interacting with a nonequilibrium quantum field was obtained. 

We show that for such a system application of the stochastic limit allows to obtain an interesting 
effect of inversion of population: for a special choice of the state of reservoir the system relaxes to 
stationary state where the population of the level of the system with higher energy will be larger 
^ ! than the population of the level with lower energy. 

We obtain an equation for the number of photons emitted and absorbed by the system. 
We investigate the examples of 2-level and 3-level systems. We show that for a 2-level atom 
the emission in the stationary regime equals to absorption. For a 3-level system emission and 
absorption of radiation are controlled by the state of the field. We find that for a 3-level system 
in a stationary nonequilibrium state two regimes are possible: the emission and the absorption 
regime. In the emission regime the total number of quanta in the system increases, and in the 
absorption regime it decreases. These regimes are controlled by the function f3(u) in ([3]). For 
example a 3-level system with energy levels E\ < e 2 < £3 is in emission regime when 

(3(e 2 - £1) + P(ea ~ £2) > /3(e 3 - £1) 

and is in absorption regime when the opposite inequality holds. 

In the emission regime the 3-level system converts radiation with the frequency 002 into radi- 
ation with the frequencies u)\ and 003 and vice versa in the absorption regime. 

This means that this stationary state of the 3-level system gives an example of dissipative 
structure in the Prigogine sense ||. 
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For an equilibrium state of the field, i.e. when the function (3(u) is linear, the system is in 
equilibrium with radiation. 

The interaction of a quantum system with a quantum field is described by an Hamiltonian 

H = H s + H R + XHj (1) 

The system degrees of freedom are described by the system Hamiltonian H$- 
The radiation degrees of freedom are described by the Hamiltonian 

H R = J oo(k)a*(k)a(k)dk (2) 

where a(k) is a bosonic field with a Gaussian state of the following form 

(a*(k)a(k')) = N(k)5(k - k') = -^—^ {k - k') (3) 

and /3(u(k)) is a (non necessarily linear) function. 

The interaction Hamiltonian Hj is defined as follows 

Hj = J ~gJk)a{k)D*dk + h.c. (4) 

We investigate the dynamics of this system in the stochastic limit, cf. in the regime of 
weak coupling (A — > 0) and large times. This regime is given by time rescaling to t *■ p-. This 
rescaling and the interaction (Q) lead naturally to introduce the rescaled quantum fields 

1 _^(a,(fc)- w ) a ^ ) (5) 
A 

where u are the Bohr frequencies (differencies of eigenvalues of the system Hamiltonian Hs)- 

By the stochastic golden rule || the rescaled field (|5|) in the stochastic limit becomes a quantum 
white noise b w (t, k), or master field satisfying the commutation relations 

[b u (t, k), k')} = 5^27:5(t - t')5(u(k) - u)5(k - k') (6) 

and with the mean zero gauge invariant Gaussian state with correlations: 

(&*(*, k)b w! (t', k')) = 6 UJU ,2irS(t - t')S(u(k) - u)8{k - k')N(k) (7) 

(b u (t, k)bl(t', k')) = 5 u ^2ir5(t - t')5(uj(k) - u)6(k - k')(N(k) + 1) (8) 

White noises, corresponding to different frequencies u, are independent. 

The Schrodinger equation becomes a white noise Hamiltonian equation, cf . , @ which when 
put in normal order is equivalent to the quantum stochastic differential equation (QSDE) 

dU t = (-idH(t) - Gdt)U t ■ t > (9) 

with initial condition Uq = 1 and where 

(i) h(t) is the white noise Hamiltonian and dH(t), called the martingale term, is the stochastic 
differential: 

rt+dt 

dH(t) = / h(s)ds = £ (E* (D) dB w (t) + (D) dB^t)) (10) 
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driven by the quantum Brownian motions 

rt+dt r rt+dt 



rt+at r rt+at 

dB w (t):= dkg(k)b u (T,k)dT=: b UJ (T,g)dr (11) 

Jt J Jt 

ii) The operator G, called the drift, is given by 



G = E {{J9\9)ZK ( D ) ( D ) + (9\9tE u (D) E* u (£>)) (12) 
where the explicit form of the constants {g\g)^ 7 called the generalized susceptivities, is: 



-i(N(k) + l) 



(g\g)Z = dk \ g (k)\ 2 I w ' = (13) 

io(k) — uj — tO 

7T / dk\g(k)\ 2 (N(k) + l)5(tv(k) - tv) - tP.P. ( dk\g(k)\ 2 ^ N ^ + V) 

J Loyk) — to 



< </!</); jdk\g(k)\ 2 ~ tN{k) = (14) 
io(k) — uj — tO 

7r / dk\g(k)\ 2 N(k)5(to(k)-io)-iP.P. [dk\g{k)\ 2 N ^ 



Ld(k) — LU 

In the present paper we consider a generic quantum system, for which for each Bohr frequency 
uj there exist a unique pair of eigenstates \lj) and (2^) corresponding to the two energy levels, 
61^,62^, so that 

In this case 

E U (D) = (1 W |D|2 W )|1 W )<2 W | (15) 
E U (D)E* W (D) = \(1 U \D\2 U )\ 2 \1 U )(1 U \ 
E^E^D) = \(1 UJ \D\2 UJ )\ 2 \2 UJ )(2 UJ \ 

We consider a dispersion uj(k) which is > and, moreover, we suppose that the Lebesgue 
measure of the set {k : uj(k) = 0} equals to zero. This implies that the real part of generalized 
susceptivities Re (g\g)t is non-negative and can be non-zero only for to > 0. 

We will also use the notation (g\g)fj for (g\g)t ii uj = Ei — Ej. 

In the present paper we investigate the non-eqilibrium stationary states for the master equation 
satisfied by the diagonal part of the density matrix of a generic quantum system. This equation 
was obtained in HI and is 



j t p{*,t) = £ (p(a',t)2Re (g\g)-, a - p(a, t)2 Re (g\g)$ a )\{a, Da')\ 2 + 

a':e (y i>e a 

+ £ (p(a\t)2Re(g\g)l,-p(a,t)2Re(g\g)- a ,)\(a',Da)\ 2 (16) 
where p(u,t) = p(a,a,t) and \a) are eigenvectors of the system Hamiltonian H$. 
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If the system has a finite number of energy levels, there exists a stationary state for the evolution 
driven by the above master equation and, if the state of the reservoir is non-equilibrium, then the 
stationary state does not satisfy the detailed balance condition for the master equation. 

The diagonal and the off-diagonal terms of the reduced density matrix evolve separately. The 
off-diagonal part of the density matrix evolves independently and vanishes exponentially, cf. |J. 
This corresponds to the collapse of a quantum state to a classical mixed state, described by the 
diagonal part of the density matrix. The diagonal part p(er, a, t) may be considered as a classical 
distribution function and equation (|T^) may be considered as a kinetic equation it. 

The structure of the present paper is as follows. 

In section 2 we describe the stationary state for a 3-level atom interacting with radiation found 
in §. 

In section 3 we investigate the properties of this stationary and find that the quotient of the 
populations of energy levels in the considered nonequilibrium stationary state does not satisfy the 
Einstein emission/absorption relation. We find that in particular case the quotient of populations 
will satisfy a new relation that we call the Double Einstein relation. 

In section 4 we use the form of this state to describe the inversion of population in our 3-level 
system. 

In section 5 we derive a master equation for the density of photons. 

In section 6 we use this equation to investigate emission and absorption of radiation by the 
system in the non-equilibrium stationary state, obtained in J3l. 



2 Stationary state for 3— level system 

Consider a 3-level system with energy levels |1), |2), |3) with energies e± < 62 < £3 and Bohr 
frequencies 

LUl = £2 — £l, ^2 = £3 — £l, <^3 — £3 — £2 

In the work || we found that the dynamics generated by the master equation (equation for 
reduced density matrix of the system) describes relaxation of the system to a stationary state. In 
this state the off-diagonal (in the basis of egenvectors of the system Hamiltonian Hg) elements of 
the reduced density matrix are equal to zero (this describes the effect of decoherence of a system 
coupled to reservoir) and the diagonal elements of the reduced density matrix for 3-level system 
take the form: 

+Kli»>ri<203)l' 1 _ / ^ M)e ^ 1 (17) 

2|/ono\|2 J M IM / 1Q x 



+ |(LD3)| 2 |(2 J D3>| 



/3(w 2 ) _ 1 1 _ e -/%a) 



e 



P3 = \(lD2)mD3)\ \ 4g0 . + |(1D2)H(2D3)| 2 ' M ' M 
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I(u) = / \g(k)\ 5(u(k) - u)dk 



where pi = (i\p\i) and 



3 The Double Einstein formula 

Consider the quotient 

Re (g\g)z = N(u) + 1 

Re (<?!<?)+ N(u) 1 ] 

Recalling that N(uj) is the density of field quanta (photon, phonons,. . . ) at frequency u, and 
comparing formula (|2(]) with the well known formula of radiation theory (the Einstein formula) 

W ■ ■ n + 1 

1 5 omission ' ^21) 



^absorption 



n,. 



giving the quotient of the probability of emission and absorption of a light quantum by an atom 
(cf. ||, Chap. V, paragraph 17, formula (18)), we gain some physical intuition of the meaning 
of the generalized susceptivities. In fact the quotient ( |2T|) . . .is just that which is necessary to 
preserve the correct thermal equilibrium of the radiation with the gas . . . (0, p. 180). 

In the stochastic limit approach this statement can be proved using the master equation (O) 
for the diagonal part of density matrix. One can prove that, if the state of reservoir is equilibrium, 
then the dynamics generated by the master equation describes the relaxation to equilibrium state 
of the system satisfying the detailed balance condition for (fLED i.e.: 



Pa, _ Re (g\g)z _ N(u>) + 1 
Pa Re {g\g)Z N(u) 

For equilibrium state the quotient of populations of the two levels with energy difference u is equal 
to the Einstein emission-absorption quotient for quanta with energy to. This suggests that the 



quotients (|20"D may play a similar role for some stationary non equilibrium states. 



Let us give an example of such a state for which we will get a generalization of condition (|22|) . 
Consider the state (O), (|18D, ([19]). For simplicity we consider the case when the matrix element 
(1-D2) is negligible (direct transitions between levels 1 and 2 are prohibited). In this case 

P2 = Re (g\g)+ 2 Re (g\g)Z 3 = N{u 2 ) N{u 3 ) + 1 = 0M j, iU3) , , 

Pl Re (g\g)- 2 Re (jg\g)+ N(lu 2 ) + 1 N(u 3 ) ' 1 ) 

Comparing with (j22D and the Einstein emission-absorption relation we will call this formula the 
Double Einstein formula. 

Relation (^) for the considered system is natural, since direct transitions from level 2 to level 
1 are prohibited ((1D2) is negligible). In this case to jump from level 2 to level 1 the system have 
to make two consequent jumps: from level 2 to level 3 and then from level 3 to level 1. Therefore 
it is natural to represent (Bl) in the following form: 



Pi 



P2 W. 



^absorption 



emission 



W ■ ■ 
yv emission 



2-3 ^absorption 



N(u 2 ) + 1 N(u 3 ) 
i_3 ~ N(u 3 ) + 1 
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Note that this formula is true for a special choice of the system ((1D2) = 0). Moreover, for a 
Gibbs state, (3 is linear and (^) coincides with the Einstein relation with N ^J^ at the RHS. 

4 Inverse population state 

Let us consider the following question: for a Gibbs distribution we have pi> p 2 > Pz'- the number 
of particles at an energy level decreases with increasing of energy Can we find such a stationary 
state where at least one pair of levels has inversed order: the number of particles increases with 
increasing energy? Such kind of states are important in quantum optics (laser theory). 

Let us apply the method of the previous section to construct a stationary state where p 2 > p\ 
(population of level 2 is larger than population of level 1). 

We will consider the same system considered in the previous section. In particular, we take 

(1D2) = 

Then from (^) we have p 2 > pi if and only if 

e -/J(w2)g/9(w8) > l 

This inequality is equivalent to: 

/3(u 3 )>P(u 1 + u 3 ) (24) 

This means that the local temperature function is non-monotonic and can decrease with an 
increase of energy. Let us note that the quotient pij p\ in the considered approximation (1D2) = 
(metastable level 2) does not depend on (1-D3) and (2D3). We found that for non-monotonic 
temperature functions we can have the inverse population effect: population of level with higher 
energy is larger than population of level with lower energy. In the theory of lasers the inverse 
population effect sometimes is discussed as an effect of negative temperature 0. Indeed if we 
suppose that the state of the field is equilibrium and therefore the local temperature function is 
linear (3{uo) = (3uj then (|24]) takes the form 

put < 

for uj\ > 0. 

In our approach we found the inverse population effect without introduction of any negative 
temperature. This effect follows from the fact that the reservoir is highly non-equilibrium and 
the temperature function can decrease with energy. 

5 Master equation for the number operator 

Consider the number operator n(k) = a*(k)a(k). This operator has constant free evolution 

e itHo n(k)e- ltHo = n(k) 

and therefore does not change in the stochastic limit. The relation with the master field is as 
follows 

[b u (t,k),n(k')] = lim^e"^^^), «(*')] = 
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= lim \e~^ {uj{k) ^ ] a{k)5{k - k') = b u {t, k)8{k - k') 
A— »o A 

This means that the number operator extends the quantum noise algebra. 

Let us find the master equation for the number operator. Since number operator does not 
commute with the noises we can not directly apply the master equation from ||. Applying the 
QSDE for the evolution operator we get 

d{U*n(k)U t ) = J2(U; (K (D) dB„{t)n{k)E u (D) 



+Eu (D) dB*(t)n(k)E* (D) dB u (t) - dtn(k)2Re GjU t ) = 
= IKfe ( D ) E » i D ) [dB U] (t),n(k)]dK(t)+ 
(D) E: (D) [dBl(t)Mk)}dB.(t)y t ) = 

■ dtJ2(U: (K P) K (D) \g(k)\ 2 2Tr5(u;(k) - oo)(N(k) + 1)- 
-E u (D) El (D) \g(k)\ 2 27r5(uj(k) - u)N(k))u t ) 



Denoting 

(X) t = (U*XU t ) 



±(n(k)) t = 2^{u{k)~u)(\g{k)\\N{k) + l)(E*(D)E ul (D)) t 



we arrive to master equation 

d_ 
dt 

-\g(k)\ 2 N(k)(E„(D)E:(D)) t ^) (25) 

This is a completely general master equation for the number operator n(k). Therefore to find 
(n(k)) t it is sufficient to determine (E*(D)E L ,(D)) t and (E u (D)E*(D)) t . 
Since we consider a generic system by flip]) equation ( |25| ) takes the form 

j t (n(k)) t = 2rr\g(k)\ 2 j:^(k) - W )|(l w |D|2 w >| 2 ((iV(fc) + l)p 2u - N(k) Pl J (26) 

6 Interaction of atom with radiation in the stationary 
state 

For a 2-level atom there is only one term in the a;-summation in (p6|) and the stationary state is 
given by (up to normalization) 

pi = Re (g\g)~, p 2 = Re{g\g)+ 
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and therefore the quotient satisfies the Einstein relation: 

e^ w) (27) 



Pl _ N(u) + 1 _ Mw) 



P2 N(u) 

If N(k) = N(u>(k)) for stationary state for two-level system (|26|) takes the form 
d 



dt (n(k)) t = 2n5(uj(k) - u)\g(k)\'\(l u \D\2j)\ 2 \N(u)(N(u) + 1) - (N(u) + l)N(u)j = (28) 

that means that in the stationary state the atom is in equilibrium with the field. Vanishing of 
(|28| ) is equivalent to the Einstein relation. Let us note that we did not assume that the state of 
the field is equilibrium. 

For a 3-level system the master equation for the number of photons (J2SI) takes the form 



j t (n(k)) t = 2n\g(k)\ 2 (5(uj(k) - Wl )|(l|D|2>| 2 ((JV(wi) + l)p 2 - N(uj 1 ) Pi ) + 

+6(uj(k) - cu 2 )|(l| J D|3)| 2 ((N(lo 2 ) + l)p 3 - N(u 2 ) Pl ) + 

+5(u(k) - ^ 3 )|(2| J D|3)| 2 ((N(u 3 ) + l)p 3 - A r (^ 3 )p 2 )) (29) 

This equation describes the balance of radiation with the 3-level atom. Consider now the station- 
ary state of the atom. When the stationary state is equilibrium then in the RHS of ( p9]) each term 
vanishes separately. This implies that in an equilibrium state the system is in detailed equilibrium 
with the radiation: emission equals to absorption for every frequency. 

When the stationary state is non-equilibrium then the Einstein relation for populations of the 
levels will be not satisfied. This implies that each term in (29~) will not vanish and the atom will 
emit and absorb radiation. 

In this case (IH, ©, flU) imply 



|(l|D|2)| 2 ((iV(^ 1 ) + l)p 2 - iV(^)pi) = 

p P(un)-J3{u>2)+i3(u> 3 ) _ i 

= Kl|D|2)| 2 |(lD3)| 2 |(2D3)| 2 /(. 2 )/(. 3 ) (e ^ i) _ _ e _ PM){ePM _ 1} (30) 
|(1|£>|3)| 2 ((JV(o; 2 ) + l)p 3 - N(u 2 ) Pl ) = 

1 _ e -/3M+/?(wi)+/3(a>3) 

= |(l|D|3)| 2 |(lD2)| 2 |(2D3)| 2 /K)/(^3) (1 _ e _^ 2))(e ^ l) _ 1){ePM _ 1} (31) 
K2IDI3)! 2 ((JV(w 3 ) + l)p 3 - N(u 3 ) P2 ) = 

= |(2|D|3)| 2 |(lD2)| 2 |(lD3)| 2 /(^)/(^) (e ^ 3) _ 1)(e , K ) _ _ e -^ 2)) (32) 

Note that (|30]), fl3T|), ( |32|) contain the combination (3(u> 3 ) + /3(u>i) — (3{uj 2 ) that vanishes for an 
equilibrium state of the reservoir. 
Equation ( ^9] ) takes the form 

-(n{k))t = 2 7 r|^)| 2 |(l| J D|2)| 2 |(l J D3)| 2 |(2 J D3)| 2 



dV y " l " lx n (e/ 3 (-i)-l)(l-e^(^))( e /3(-3)_l) 



IM I{Lo 3 )5{u{k) - wi) - 1(^)1(^)5^) - u 2 ) + /(a;i)/(a; 2 )(J(w(A;) - cu 3 )J (33) 

Equation ( p3[ ) shows that the 3-level system in non-equilibrium stationary state converts radiation 
with the energy u 2 = 0J\ + u> 3 into radiation with energies oj\ and u 3 if 

PM + PM > /%i + w 3 ) 

and vice versa in the case of the opposite inequality. 
Integrating fl3"3]) over we get 



d 
dt 



E(u(k)){n(k)) t dk = 27r|(l|D|2)| 2 |(l J D3)| 2 |(2 J D3)| 2 



p P(un)-P(c*>)+P(u>a) _ i 
IMIMIM (EM - EM + EM) (e ^ 1 )-i)(i- e -^))(e^3)-l) 

When E{oj) is the dispersion of the field 

E(u) = u 

then 

EM - EM + EM = o 
and (|34]) implies the conservation of energy 



u(k)(n(k)) t dk = 



For the time derivative of number operator we get 



j t {n(k)) t dk = 27r|(l|D|2>| 2 |(lD3>| 2 |(2 J D3>| 

e P(ui)-P(U2)+P{u)3) _ i 

iMH^iM- 



(ePW - 1)(1 - e-^))(e^( u s) - 1) 
Using that cu 2 = c^i + we obtain that if 

PM + PM < P(^ + ^3) (35) 

then the derivative is negative and the system absorbs the radiation (total number of absorbed 
photons is larger than the total number of emitted photons). In this case (j33|) implies that the 
system absorbs photons with frequencies ui\ and uj 3 and emits photons with frequency uj 2 . 
If 

PM + PM > P(ui + 003) (36) 

then the derivative is positive and the system emits the radiation (the total number of absorbed 
photons is smaller than the total number of emitted photons). In this case fl33|) implies that the 
system emits photons with frequencies ui\ and u 3 and absorbs photons with frequency u 2 . For 
instance in the case of inverse population (EM) the condition of emission regime (59) is satisfied. 
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These regimes of emission and absorption are controlled only by the difference (3(ui) + (3{uj$) — 
/3(a>i + U) 3 ). 
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